1/ f a noise in spectral fluctuations of quantum systems 
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The power law l/f a in the power spectrum characterizes the fluctuating observables of many 
complex natural systems. Considering the energy levels of a quantum system as a discrete time 
series where the energy plays the role of time, the level fluctuations can be characterized by the 
power spectrum. Using a family of quantum billiards, we analyze the order to chaos transition in 
terms of this power spectrum. A power law l/f a is found at all the transition stages, and it is 
shown that the exponent a is related to the chaotic component of the classical phase space of the 
quantum system. 
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Of the different features which characterize complex 
physical systems, perhaps the most ubiquitous, interest- 
ing and puzzling is the presence of 1//" noise |l| in fluc- 
tuating physical variables, i.e. the Fourier power spec- 
trum S(f) behaves as l/f a in terms of the frequency /. 
This kind of noise has been detected in condensed matter 
systems, traffic engineering, DNA sequence, quasar emis- 
sions, river discharge, human behavior, heartbeat and dy- 
namic images, among many others. Despite this ubiquity, 
there is no universal explanation about this phenomenon. 
It does not arise as a consequence of particular physical 
interactions, but it is a generic manifestation of complex 
systems. 

Recently, it was conjectured that the energy spectra of 
chaotic quantum systems are characterized by 1/f noise 
. The original idea was that the sequence of discrete 
energy levels in a quantum system can be considered as 
a discrete time series, where the energy plays the role 
of time. In that case, the energy level fluctuations can 
be studied using traditional methods of time series anal- 
ysis, like the study of the power spectrum. When the 
idea was applied to typical chaotic quantum systems, the 
power spectrum showed a very accurate 1/f behavior Q • 
Hence chaotic quantum systems can be added to the long 
list of complex natural systems which exhibit 1 // noise. 
However, this new point of view also raises new ques- 
tions. Is this a consequence of the universal behavior of 
fluctuations in chaotic quantum systems? What happens 
in quantum systems which are neither fully chaotic nor 
fully regular? In this paper we try to find some answers 
using a quantum billiard to study the power spectrum in 
the order-to-chaos transition. As shown below, the ubiq- 
uitous 1// Q noise appears at all the transition stages, 
with the exponent smoothly decreasing from a = 2 in 
a regular system to a = 1 in a chaotic system. This is 
quite a remarkable result indeed, since it contradicts the 



predictions of the strict semiclassical limit |i| • 

The concept of quantum chaos, or wave chaos in more 
general terms 0, has no unique precise definition as yet, 
but definitely can be described as quantum or wave like 
signatures of classical chaos. It is well known that there is 
a relationship between the energy level fluctuation prop- 
erties of a quantum system and the dynamics of its clas- 
sical limit. Classically integrable systems give rise to 
uncorrelated adjacent energy levels, which are well de- 
scribed by Poisson statistics 0- In contrast, spectral 
fluctuations of a quantum system whose classical limit is 
fully chaotic (ergodic) show a strong repulsion between 
energy levels and follow the predictions of random ma- 
trix theory (RMT) [10. In practice, quantum systems 
without classical limit are assumed to be chaotic when 
their fluctuations coincide with RMT predictions. 

The essential feature of chaotic energy spectra is the 
existence of level repulsion and correlations (leading to 
strong spectral rigidity), i. e. the spacing of two ad- 
jacent levels is unlikely to deviate much from the mean 
spacing. This property is similar to the antipersistence 
characteristic of some time series Q. Antipersistence, 
with different intensity degrees, appears in time series 
with l/f a noise, with I < a < 2. Could the analogue 
level repulsion feature be also associated to I / f a noise? 

To study the spectral fluctuations of quantum systems 
we follow the method introduced in Q ■ We use the statis- 
tic S n denned by 

n 

Sn = ^ ( Si ~ ( S )) = 6,1+1 ~ 61 ~ U > W 
i=l 

where are the unfolded energy levels, Sj = e.; + i — e,, 
and < s >= 1 is the average value of Si. Thus 5 n repre- 
sents the fluctuation of the nth excited state. Formally 
S n is similar to a time series where the level order index 
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n plays the role of a discrete time. Therefore the sta- 
tistical behavior of level fluctuations can be investigated 
studying the power spectrum S(k) of the signal, given by 



S(k) = 
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where M is the size of the series and / = 2nk/M plays 
the role of a frequency. 

To investigate the behavior of S(k) in the mixed regime 
between integrability and chaos, we analyze it in the Rob- 
nik billiard |g. Quantum billiards are considered as a 
paradigm in quantum chaos. They have a discrete spec- 
trum with an infinite number of eigenvalues, and there- 
fore it is possible to reach high statistical precision by 
computing a large number of them. Furthermore, they 
can also be studied experimentally 0, . 

The boundary of the Robnik billiard is defined as the 
set of points w in the complex plane C which satisfy 
the equation w = z + Xz 2 , where \z\ = 1 and A is the 
deformation parameter. It has been shown || that this 
billiard exhibits a smooth transition from the integrablc 
case (A = 0) to an almost chaotic case (1/4 < A < 1/2). 
In order to obtain a smooth analytic boundary, A must lie 
in the interval [0, 1/2). The Robnik billiard is one of the 
be st sy stems to investigate the order-to-chaos transition 
|sL lld lTl) . Compared to other quantum billiards, it has 
the advantage that there are no bouncing ball orbits. For 
small values of A the billiard is a typical KAM system, 
whereas for larger values of A only one chaotic region 
dominates the phase space with only few stability islands 
covered with invariant tori. The total area in the bounce 
map (Poincare surface of section) of these invariant tori 
decreases monotonically with A and becomes negligible 
when the shape of the billiard becomes non-convex, for 
A > 1/4. For A = 1/2 it has been shown rigorously by 
Markarian that the billiard is ergodic. 

The quantum energy levels E n of the Robnik bil- 
liard are numerically calculated by solving the station- 
ary Schrodinger equation of a free particle whose wave 
function ip{w) is zero at the boundary of the billiard. 
The billiard has reflection symmetry with respect to the 
real axis, so there are two types of states: those with 
even parity ip(w) — ip(w*) and those with odd parity 
ip(w) — —ij){w*)\ odd and even parity states must be 
treated separately 0, 0] . For each symmetry, our cal- 
culation uses approximately 80,000 basis states, giving 
good eigenvalues for about 65,000 levels for A = and 
about 30,000 levels for A = 0.5. 

Fig. ^ shows the energy level fluctuations of the Rob- 
nik billiard given by S n . It illustrates the effect of level 
repulsion in the order-to-chaos transition, and its rela- 
tionship with the antipersistence of S n considered as a 
time series. For the regular system (A = 0), the levels 
are uncorrclatcd and therefore S n is neither persistent 
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FIG. 1: Plot of the statistic S n for a set of 256 consecutive 
energy levels of odd parity in the Robnik billiard, for several 
values of the deformation parameter A. 



nor antipersistcnt. As A increases the system becomes 
more chaotic, and S n looks like a typical antipersistent 
series in the almost chaotic region for A > 1/4. 
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FIG. 2: Average power spectrum (S(k)) of the statistic S n for 
the odd parity energy levels corresponding to the shapes of 
the Robnik billiard inserted in the figures. The four values of 
the deformation parameter A are the same as in Fig. Q The 
solid (red) line is the best fit to the power law l/k a . 

Fig. El shows the power spectrum of S n obtained from 
the energy levels E n of the Robnik billiard for several 
values of A. The contour of the billiard is shown as an 
inset in each panel. We calculate an ensemble average of 
S(k) in order to reduce statistical fluctuations and clarify 
its main trend. The average (S(k)) is calculated with 25 
sets of 256 consecutive high energy levels of odd parity. 
It is clearly seen that for each value of A it follows the 
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scaling law 
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where a depends on A. In fact the fit of (S(k)) to the 
power law l/k a is excellent. In all cases the error in 
the linear regression is less than 3%. For A = (inte- 
grable case) the exponent is a = 1.98, as expected for 
uncorrelated energy levels. As A increases the exponent 
a decreases and becomes a ~ 1 for A ~ 1/2. Thus, a 
may serve as a measure of the chaoticity of the system, 
since it changes from a = 2 for regular systems to a = 1 
for chaotic ones. 

It is worth to compare this result with the behavior 
of more conventional statistics, like the nearest neighbor 
spacing distribution P(s), which measures short range 
correlations, and the Dyson A^(L) statistic appropriate 
for correlations of length L [l^. Fig. ^displays P(s) for 
several values of A. At A = the histogram follows the 
predicted curve for regular systems (Poisson limit). For 
A = 0.15, P(s) deviates from Poisson toward the RMT 
limit. As we shall see below, this behavior of P(s) reflects 
that the underlying classical dynamics is neither regular 
nor ergodic (chaotic). Finally, for A = 0.25 and A = 0.4 
the system exhibits short range correlations characteris- 
tic of chaotic systems (RMT limit). Fig. 0] shows the 
spectral average (A3(L)) for energy intervals of length 
L ranging from L = 2 to L = 50, for several values of 
A. The spectral average is calculated using 25 sets of 
L consecutive high energy levels to avoid, as far as pos- 
sible, the influence of short periodic orbits. The evolu- 
tion of this statistic with A is analogous to that of P(s). 
When A = 0, (A(L)} falls near the Poisson prediction for 
regular systems, and for A = 0.25 and 0.4 it is almost 
indistinguishable from the RMT prediction for chaotic 
systems. Therefore, P(s) and (Aa(L)) have a smooth be- 
havior in terms of A. Nevertheless, they move faster than 
S n toward the RMT limit as A increases. For instance, 
both P(s) and (A 3 (L)) coincide with RMT predictions 
for A = 0.25, while S n still points to an intermediate 
regime between regularity and chaos. 

Let us now compare the evolution of the parameter a 
as a function of A with the fraction pf of regular classi- 
cal trajectories in the phase space, and with the Brody 
parameter u> |l4j . This is an ad hoc parameter, without 
any known physical meaning, which quantifies to some 
extent the chaoticity of the system. It nevertheless cap- 
tures well the important feature of fractional power law 
level repulsion, that is the behavior of P(s) at small s 
fiot llH . For w=0we get the Poisson distribution of reg- 
ular systems, and at the other extreme, u) = 1, we obtain 
the Wigner distribution predicted by RMT for chaotic 
systems. Fig. shows the behavior of a, u and pf. 
There is a clear correlation among these three variables, 
although the transition is smoother for oj and especially 
for a than for pf. However, while the fraction of regu- 
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FIG. 3: Nearest neighbor spacing distributions for the spectra 
of four different shapes of the Robnik billiard. The solid line 
is the Poisson distribution and the dotted line corresponds to 
the Wigner distribution predicted by RMT. 
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FIG. 4: Spectral average (Aa(L)) values calculated numeri- 
cally for the Robnik billiard using the same deformation pa- 
rameters as in previous figures. The results of the Poisson 
statistics (solid line) and RMT (dotted line) are also dis- 
played. 



lar classical trajectories is almost zero near A = 0.25, the 
power spectrum, and to a lesser extent the P(s) distribu- 
tion, indicate an intermediate situation between regular 
and chaotic motion. This clearly shows that w and a are 
not only functions of pf, but depend on finer details of 
the underlying classical mechanics as well. 

It is well known that in the strict scmiclassical limit the 
quantum cigenstates of a quantum system with generic 
(mixed) classical dynamics can be classified as regular 
and irregular following the original proposition of Perci- 
val . This has been further developed and raised to a 
Principle of Uniform Semiclassical Condensation (PUSC) 
of Wigner functions of the eigenstates From this 

it follows that in the strict semiclassical limit the reg- 
ular and irregular level sequences are statistically inde- 
pendent, but for themselves have Poisson or RMT level 
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provides an intrinsic quantitative measure of the regular 
and chaotic dynamical features. 

The origin of the universal power law behavior S(k) ~ 
l/k a is now understood in the intcgrable case (a = 2) 
and in the fully chaotic case (a = 1) on the basis of RMT 
[l5j and semiclassical periodic orbit theory Q . The origin 
of the 1//" power law in the mixed regime still remains 
as an important open problem. 

This work is supported in part by Spanish Govern- 
ment grants BFM2003-04147-C02 and FTN2003-08337- 
C04-04. This work is also supported by the Ministry of 
Education, Science and Sports of the Republic of Slove- 
nia. 



FIG. 5: Behavior of the power spectrum exponent a, the 
Brody parameter lo, and the fraction pf of regular classical 
trajectories of the Robnik billiard as functions of the defor- 
mation parameter A. 

statistics, respectively. This theory has been excellently 
confirmed in hard numerical calculations 01 ■ Neverthe- 
less, if the system is not deep enough in the semiclassical 
regime one can see substantial deviations from such a be- 
havior, manifested in the fractional power law level repul- 
sion [l(j,lllj in P(s) at small s. A similar recent analysis 
has demonstrated that in the strict semiclassical limit 
we should not expect a power law behavior for the power 
spectrum but something more complicated. Therefore it 
is quite an unexpected result of the present paper that 
the power spectrum S(k) is a power law at all k in mixed 
systems at low energies. Indeed, when going sufficiently 
deep into the semiclassical regime the theory of reference 
should be expected and confirmed. 

In conclusion, the analogy between quantum energy 
spectra and time series opens a new and fruitful perspec- 
tive on the universal properties of quantum level fluctu- 
ations. The 5 n function gives the level fluctuations, and 
its power spectrum S(k) is an intrinsic characteristic of 
the quantum system. The important point is that it ex- 
hibits a power law behavior, similar to the well known 
1// Q noise found in many complex systems. 

In the order to chaos transition, the chaoticity of a 
quantum system is usually qualitatively assessed by how 
close to Poisson or RMT its fluctuation properties are. 
In the present power spectrum approach, the exponent 
changes smoothly from a = 2 for a regular system to 
a = 1 for a chaotic system. Contrary to the Dyson A^(L) 
statistic, that must be plotted for different values of L, 
the exponent a quantifies the chaoticity of the system in 
a single parameter. Moreover, a has a physical mean- 
ing. It is a natural measure of the fluctuation properties 
of a quantum system through the power spectrum, and 
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